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Abstract. The fusion-fission process accompanied by neutron emission is studied in 
the superheavy-mass region on the basis of the fluctuation-dissipation model combined 
with a statistical model. The calculation of the trajectory or the shape evolution in 
the deformation space of the nucleus with neutron emission is performed. Each process 
(quasi-fission, fusion-fission, and deep quasi-fission processes) has a characteristic 
travelling time from the point of contact of colliding nuclei to the scission point. 
These dynamical aspects of the whole process are discussed in terms of the pre-scission 
neutron multiplicity, which depends on the time spent on each process. We have 
presented the details of the characteristics of our model calculation in the reactions 
^^Ca-|-^°^Pb and ^*Ca-|-^'*'*Pu, and shown how the structure of the distribution of 
pre-scission neutron multiplicity depends on the incident energy. 
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1. Introduction 

Recently, spectacular news on the finding of superheavy elements has been successively 
reported from FLNR, GSI and RIKEN [H H El H E]. In theoretical approaches, 
the clarification of the reaction mechanism in the superheavy mass region and the 
prediction of the favorable conditions (projectile-target combination and beam energy 
values) to synthesize new elements have advanced step by step in the last few years 
[3 ISl El Uni HH nil Uni- However, we are still seeking suitable models for calculating 
the fusion and the survival probabilities, and investigating the unknown parameters 
involved in them. The fusion-fission mechanism in the super heavy-mass region must 
be elucidated; in particular, a more accurate estimation of fusion probability should be 
estabhshed. 

Generally, on the base of analysis of the experimental mass, kinetic energy, and 
angular distributions of the reaction fragments, the whole reaction process in heavy 
nucleus collisions is classified into the fusion-fission, quasi-fission and deep inelastic 
collision processes and so on. Generally, mass symmetric fission fragments have been 
considered to originate from a compound nucleus. Therefore, in the experiment, the 
fusion-fission cross section has been derived by counting mass symmetric fission events 

However, in our previous work we showed that it is ambiguous to identify mass 
symmetric fission events as the fusion-fission process in the superheavy-mass region. 
Using the Langevin equation, which enables to calculate the time development of a 
trajectory of a colliding system in a shape parameter space, we attempted to distinguish 
different dynamical paths by analyzing the trajectories and mass distribution of fission 
fragments. We identified the fusion-fission (FF) path, the quasi-fission (QF) path and 
the deep inelastic collision (DIG) path. In addition, in the reaction ^^Ga+^^^Pu, we 
found a new path, called the deep quasi-fission (DQF) path, whose trajectory does not 
reach the spherical compound nucleus but goes to the mass symmetric fission area. 
The conclusion is that the analysis of only the mass distribution of fission fragments is 
insufficient to identify the reaction process described between the FF and DQF paths. 
To estimate the fusion probability accurately, the precise identification of the FF path 
is very important. 

On the basis of experimental results, it has also been suggested that the reaction 
process may be classified by the measurement of the pre-scission neutron multiplicity, 
which correlates with the mass distribution of fission fragments [T^ [TB| IT7j . This means 
that each process is expected to have its own characteristic reaction time, which is 
related to the pre-scission neutron multiplicity. 

Hence, to identify the reaction process more precisely, we undertake to extend our 
model discussed in reference [12], by taking into account the effect of neutron emission 
along the dynamical path. We introduce the effect of neutron emission into the three- 
dimensional Langevin calculation code, or we combine the fluctuation and dissipation 
dynamical model with a statistical model. The preliminary analysis along this line has 
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been done in the reaction system ^^Ni+^^'^Pb at a incident energy corresponding to the 
excitation energy of the compound nucleus E* = 185.9 MeV, which is a high incident 
energy case [TH] . 

In this study, we have apphed our model to the investigation of the whole 
dynamical process in the reactions ^^Ca+^°^Pb and ^^Ca+^'^^Pu mainly at E* = 50 
MeV, that is to say, a low incident energy case. These reactions are chosen in order 
to precisely investigate the fusion process in the hot fusion reaction that has recently 
been performed at FLNR ^H]- Also we discuss how the characteristic structure of the 
neutron multiplicity depends on the incident energy. 

In section 2, we explain our framework for the model used here. We discuss the 
pre-scission neutron multiplicity with respect to fission fragments in section 3. The 
mechanism of the fusion-fission process with neutron emission is investigated. In section 
4, we present a summary and further discussion. 



Using the same procedure as described in reference J3] to investigate the whole 
process dynamically, we use the fluctuation-dissipation model and employ the Langevin 
equation. We adopt the three-dimensional nuclear deformation space given by two- 
center parameterization fl^ 1^ . The three collective parameters involved in the 
Langevin equation are as follows: zq (distance between two potential centers), S 
(deformation of fragments) and a (mass asymmetry of the colliding nuclei); a = 
{Ai — A2)/{Ai + A2), where Ai and A2 denote the mass numbers of the target and 
the projectile, respectively. We assume that each fragment has the same deformations 
as the first approximation. In the present calculation, the neck parameter e is fixed to 
be 1.0, so as to retain the contact-like configuration more realistically for two-nucleus 
collision. The definitions of the 6 and e are mentioned in reference fl^ 121( 1^. 
The multidimensional Langevin equation is given as 



where a summation over repeated indices is assumed, qi denotes the deformation 
coordinate specified by Zq, S and a. Pi is the conjugate momentum of q^. V is the 
potential energy, and mij and are the shape- dependent collective inertia parameter 
and dissipation tensor, respectively. A hydrodynamical inertia tensor is adopted in the 
Werner- Wheeler approximation for the velocity field, and the wall-and-window one-body 
dissipation is adopted for the dissipation tensor [231 1211 • The detail is explained in 
reference |15j . 

The intrinsic energy of the composite system Eint is calculated for each trajectory 

as 



2. Model 
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where E* denotes the excitation energy of the compound nucleus, and is given by 
E* = Ecm — Q with Q and Ecm denoting the Q— value of the reaction and the incident 
energy in the center-of-mass frame, respectively. 

We take into account neutron emission in the Langevin calculation during the 
reaction process. The emission of neutrons which is calculated by the statistical 
model has been coupled to the three-dimensional Langevin equation. The proceeder 
is explained precisely in reference |18j . 

3. Reaction process and pre-scission neutron multiplicity 

Our model is applied to the reaction systems ^^Ca+^°^Pb and ^^Ca+^^^Pu with the 
incident energy corresponding to the excitation energy of the compound nucleus, E*=50 
MeV. In our previous work pH], we have shown that our model calculations for the 
mass distribution of fission fragments agree well with the experimental data and also 
the fusion-fission cross section. On the basis of the results given in reference ^3], we 
want to make a new trial to clarify the relationship between the dynamical process 
and the pre-scission neutron multiplicity, and to investigate whether the relation can 
be seen at E*=50 MeV. These reactions are chosen bearing in mind the more precise 
estimation of the fusion probability in the so-called hot fusion reaction performed at 
FLNR in experiments on superheavy element synthesis jl9j . 

In the preliminary investigation in reference ^Hl, for the purpose of checking the 
applicability of our model, we have analyzed to the recent experiment, in which the pre- 
scission neutron multiplicity correlated with the mass distribution of fission fragments 
has been measured in the reaction ^^Ni+^°^Pb at the incident energy corresponding to 
the excitation energy of compound nucleus E* = 185.9 MeV. This experiment was done 
by DeMoN group ^Hl- Our model calculation has reproduced the characteristic trend 
of the experimental distribution of pre-scission neutron multiplicity, and distinguished 
the FF and QF processes among the mass symmetric fission events. 

Also we have showed our preliminary results on this study in references |17| |2H1 
EZlEHj- We investigated the neutron multiplicity in the reaction ^^Ca+^^^Pu at E*=4:0 
MeV. The experiment of neutron multiplicity shows the two peaks near 2 and 4 neutron 
emission ^7j. Such two peaks are originated by two different fusion-fission mechanism. 
It looks that the first peak comes from QF process and the second one comes from mainly 
DQF or FF. However, the calculation could not reproduce the two peaks j2H]- Due to a 
low excitation energy, the lifetime of neutron emission is rather short. It is comparable 
to the time scale of the fluctuation of the trajectory. So we can not distinguish QF and 
DQF process clearly. Moreover the number of events of FF process is quit small at a 
low incident energy. The statistics of each event are not enough to analyze the process. 

Here, we investigate further by applying in the reaction '^^Ca+^°'^Pb and 
^^Ca+^^^Pu, more precisely, with increasing the excitation energy E* = 50 MeV. To 
discuss the dynamical process, we mainly focus on zero angular momentum case. For 
each angular momentum, the potential landscape changes but the mechanism of the 
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process essentially does not change [221 • 



3.1. Reaction ^^Ca+^^^Pb 

Using the model presented in Sect. 2, we calculate the pre-scission neutron multiplicity 
emitted from the reaction system ^^Ca+^°^Pb at E* = 50 MeV for angular momentum 
/ = 0. In this system, the FF process is dominant, that is to say, the mass symmetric 
fission events are dominant. 




Figure 1. Pre-scission neutron multiplicities given by our theoretical calculation in 
the reaction ^^Ca+^^^Ph at E*^50 MeV for / = 0. The neutron multiplicity from the 
QF, DQF and FF processes are denoted by the dashed-dotted, the dashed and the 
dotted lines, respectively. The solid line denotes the total process. 

The calculation is done in the same procedure as reference |22t IT8] . In the three- 
dimensional Langevin calculation, we prepare 1,000 trajectories for each case. At t = 0, 
each trajectory starts from the point of contact. We define the fusion box as the 
region inside the fission saddle point, {z < 0.8,6 < 0.3, |a;| < 0.3} jTH]. Samples of 
the trajectory of the FF process projected onto z — a {6 = 0) plane in the reaction 
^^Ca+^ospb at E* = 50 MeV for / = is shown in Fig. 1 in reference |2ZI. The FF 
trajectory makes up 96% of all trajectories, and 95% of the FF trajectory escapes from 
the spherical region within 5.0 x 10~^^ sec. 

Figure 1 shows the distribution of the pre-scission neutron multiplicity in this 
reaction. The multiplicities from the QF, DQF and FF processes are denoted by the 
dashed-dotted hne, dashed line and dotted line, respectively. The solid line shows the 
total multiplicity of each process. The distribution of pre-scission neutron multiphcity 
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Figure 2. The distribution of travelling time ttrav in the reaction ^^Ca+^'^^Pb at 
E* — 50 MeV for I = 0. The tfrav from the QF, DQF and FF processes are denoted 
by the dashed-dotted line, dashed line and dotted line, respectively. The solid line 
denotes the total process. 

has a peak around i/„ = 2. Since the FF process is dominant, the peak corresponds to 
the FF process. 

We define the travelhng time ttrav as a time duration during which the trajectory 
moves from the point of contact to the scission point. We present the distribution 
of the travelhng times for the QF, DQF and FF processes in Fig. 2, denoted by the 
dashed-dotted fine, dashed fine and dotted fine, respectively. The solid line shows the 
distribution of the travelling time for all processes. In the FF process, as discussed 
in reference [221; the trajectory is trapped in the pocket around the spherical region. 
Consequently, the trajectory spends a relatively long time in the pocket before it finally 
escapes from there. The travelling time of the FF process has a peak at 1.2 x 10~^° sec, 
and the average time of it is 15.3 x 10~^° sec. 

As shown in reference jlH], in the reaction system ^^Ni+^°^Pb at E* = 185.9 MeV, 
the distribution of the neutron multiplicity has two peaks which come from the QF and 
FF processes, respectively. On the other hand. Fig. 1 shows the single peak. 

For the ^^Ca+^°®Pb reaction, the FF process is dominant in this energy region, in 
consistent with the potential energy surface shown in Fig. 1 in reference j22]; in which 
we can see the steep valley in the entrance mass asymmetry and in the mass symmetric 
fission valley. This fact is consistent with the pre-scission neutron multiplicity having 
only a single peak at z/„=2, which indicates that only one process is dominant. 

In the preliminary calculation [2H1, we have compared the results with the 
experimental data at E* = 40 MeV [iZj , which shows the distribution of the pre-scission 
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neutron multiplicity in correlation with fission fragments whose mass number is grater 
than — 30 and less than + 30, where Acn denotes the mass number of the 
compound nucleus. The tendency of these results coincides with the experimental data. 

3.2. Reaction ^^Ca+^^^Pu 

In reference [T5t l22] . the whole dynamical process in the reaction ^^Ca+^^^Pu was 
classified into the QF, DQF and FF processes. At a low incident energy, the mass 
asymmetric fission events are dominant, which is the QF process ^3]. We have 
distinguished between the DQF and the FF process by analyzing whether the trajectory 
enters the fusion box or not. The analysis only of the mass distribution for fission 
fragments is insufficient to distinguish between the DQF and the FF processes |TS|. In 
the present study, we show how each process corresponds to a characteristic number of 
pre-scission neutron multiplicities at E* = 50 MeV. 

Samples of the trajectory of the QF, DQF and FF process projected onto z — a 
{S = 0) and z-6 {a = 0) plane at E* = 50 MeV for / = in the reaction *^Ca+'^^'^Pu 
are shown in Fig. 5 in reference j22I- We define the fusion box as the region inside the 
fission saddle point, {z < 0.6,5 < 0.2, |a| < 0.25} jl5j. 

3.2.1. Classification on the basis of mass distribution of the fission fragments 
Generally, it is considered that the reaction process is classified according to the mass 
distribution of fission fragments and their total kinetic energy distribution ^3] . The QF 
process is defined as the fission fragments having mass number greater than + 20 
and less than — 20. This condition is also suited to the deep inelastic collision 
process (DIG) in the present case. The FF process is defined by the condition that the 
trajectory enters the fusion box. The DQF process produces the mass symmetric fission 
fragments (greater than — 20 and less than + 20), but the system does not 
reach the spherical compound nucleus before undergoing fission. We name this method 
of classification as the "classification by fission fragments mass" (standard classification). 

As discussed in reference ^3] , we can not distinguish between the FF and the DQF 
processes only using the fission fragment mass. The difference of the both processes is 
the reaction time, that is to say, the travelling time ttrav We suppose that the neutron 
multiplicity depends on the tt^av Here, we try to investigate the correlation between 
the neutron multiplicity and each process which is separated by using the standard 
classification. 

As the result analizing 10,000 trajectories, it is found that the QF process (including 
a part of DIG), the DQF process and the FF process account for 72%, 25% and 3%, 
respectively. In Fig. 3, the neutron multiplicity from the QF, DQF and FF processes 
are denoted by the dashed-dotted line, solid line and gray line, respectively. The dashed 
line denotes the total multiplicity contributed from all processes. 

Many trajectories go to the fission area without neutron emission. The travelling 
time of such a trajectory is rather short. As can be seen in Fig. 3, it is difficult to 
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Figure 3. Pre-scission neutron multiplicity given by our theoretical calculation in 

the reaction ^^Ca+^^^'Pu at i?*=50 MeV for I = with the classification method of 
using the mass distribution of fission fragments. The neutron multiplicity from the 
QF, DQF and FF processes are denoted the dashed- dotted, the solid and the gray 
lines, respectively. The dashed line denotes the total multiplicity of each process. 



see a characteristic structure of the neutron multiphcity corresponding to the process 
classified. The total number of events for 2n emission is significant, but the QF and 
the DQF processes account for approximately equal contributions. These results show 
that there is no correlation between the neutron multiplicity and the mass distribution 
of the fission fragments in this incident energy. 

Figure 4 shows the mass distribution of fission fragments for each process. The 
fission fragments from the QF+DIC, DQF and FF processes are presented by the light 
gray, the gray and the black shading, respectively. The thin black fine denotes the total 
number of events of all processes. The fission fragments from the FF process spread 
from A — 100 to 200. Therefore, we can see that the FF process is related to both the 
mass symmetric fission and the mass asymmetric fission. In Fig. 4, two sharp peaks are 
located at A ~ 80 and A ~ 210. These events come from the trajectories that go down 
quickly along the valley at a ~ 0.42 in the potential energy surface, which is the result 
of the the shell structure of the Pb nucleus. 

The corresponding travelling times for each process arc shown in Fig. 5. The peak 
of travelling time for the QF process, which is located at ttrav ~ 7.5 x lO"^-*^ sec, is slightly 
shorter than the time at the peak of the DQF process. Until ttrav — 3.0 x lO"^'^ sec, 
the traveUing time of DQF process overlaps with that of the QF process. The average 
traveUing times of the QF, DQF and FF processes are 1.11, 1.73 and 4.83 x 10"^° sec. 
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Figure 4. Mass distribution of fission fragments in the reaction °Ca+ Pu at 

E* = 50 MeV for I = using the method of classification by fission fragment mass. 
The fission fragments from the QF+DIC, DQF and FF processes are presented by the 
fight gray, the gray and the black shading, respectively. The thin black line denotes 
the total number of events of all processes. 




Figure 5. The distribution of the travelling time ttrav of each path in the reaction 
^*Ca+^^^Pu at E* = 50 MeV for / = 0, using the method of classification by fission 
fragment mass. 
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respectively. The travelling time of the FF process is about three times as long as those 
of the QF and DQF processes. 

3.2.2. Correlation between the travelling time ttrav and the neutron multiplicity z/„ 
In the standard classification, we can not distinguish a characteristic structure of 
the neutron multiplicity corresponding to each process, in this low incident energy 
case. Here, we investigate the reason why the neutron multiplicity do not show the 
characteristic structure corresponding to the reaction process like as a high incident 
energy case [T^ . In conclusion, the cause of this failure is the use of the fission fragment 
mass for the classification of processes. We discuss this point in this subsection. 

As we discussed above, we expect that each process may have the characteristic 
reaction time. Here, at first we show the relation between the average travelling time 
(ttrav) and the mass number of the fission fragments in Fig. 6. The average travelling time 
of the mass symmetric fission process is longer than that of the mass asymmetric fission 
process. We can see clearly the correlation between the (ttrav) and the mass number 
of the fission fragments. The average travelling time for — 20 < A < + 20 is 
ranging from 1.8x10^^° to 2.1x10^^° sec. In Fig. 5, the fiuctuations of the travelling 
time of the QF and the DQF processes, which are expressed by the distributions of 
the travelling time around the most probable value for each process, are ranging in the 
width of 2.5 ~ 3.0 x 10~^° sec. These fiuctuations of the travelling time are larger than 
the difference between the average travelling time of the QF and the DQF processes. 
Therefore, it is difficult to separate the QF and DQF processes by the travelling time. 

The relation between the travelling time and each neutron multiplicity in this 
system is shown in Fig. 7. With increasing the time duration, the large neutron 
multiplicity is expected. At this incident energy, the events for In and 2n emission 
are dominant, and it corresponds to the peaks for the QF and DQF processes in Fig. 3. 
Due to the stochastic fiuctuation in the statistical model, the time taking for the events 
of In and 2n emission are overlapped significantly. We can not distinguish the reaction 
time of the both processes by the number of neutron multiplicity. 

The life time of the neutron emission and the travelling time are partly governed 
by the random number in the stochastic aspect which is contained in statistical model 
and the Langevin equation. The width of these fiuctuation in time space is overlapped 
for each process in this case. Actually, at E* = 50 MeV, the maximum number of 
evaporated neutron z/^°^ is approximately 4 ~ 5 as shown in Fig. 3. However, as shown 
in Fig. 7, peaks for the different neutron multiplicity are not well separated so as to 
identify the reaction process. As we will show later, it is coming from small z/^°^(~ 4). 

Figure 8 shows the neutron multiplicity associated with the mass number of fission 
fragments in the reaction ^^Ca+^^^Pu at E* = 50 MeV for / = 0. The two sharp peaks 
of On and In at A = 80 and 210 are caused by the trajectory which goes down along 
the valley near the Pb nucleus (or the Kr nucleus) in the potential energy surface. Such 
a trajectory quickly goes down to the fission area, so that it does not have large chance 
to emit any neutron. As can be seen in Fig. 8, there exist trajectories with a variety of 
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Figure 6. The average travelling time associated with fission fragment in the reaction 
48Ca+244pu at E* = 50 MeV for / = 0. 
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Figure 7. Time dependence of neutron multiplicity in the reaction ''*Ca+^ Pu at 
E* = 50 MeV for / = 0. The neutron multiplicity are given. 
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Figure 8. The neutron multiplicity associated with fission fragment in the reaction 
48Ca+244pu at E* = 50 McV for / = 0. 



numbers of the neutron emission in the region 90 < A < 210. The trajectories having 
different neutron multiphcities coexist in this region. With decreasing mass asymmetry, 
2n and 3n events increase because, in this region, the number of trajectories with long 
travelhng time increases as shown in Fig. 6. 

The characteristic feature observed in Fig. 8 is that each neutron multiphcity 
spreads over a wide mass region of fission fragments. It is noteworthy that the events 
with even On and In are found in the mass symmetric fission area. Therefore, it is clear 
that we can not extract the characteristic structure of neutron multiplicity by means of 
the standard classification. 



3.2.3. Classification using the additional condition on z As discussed in the previous 
subsection, the trajectories of each process spread over a wide mass region of fission 
fragments. We understand that each process has characteristic behavior. However, the 
information of the mass number of fission fragments does not clearly connected with 
the dynamical aspect of each process which we imaged. We have to look for the proper 
classification method that corresponds to our recognition of each process, instead of the 
method using fission fragment mass. Here, we propose to employ a classification method 
based on the dynamical trajectory of the process. 

As shown in Fig. 3, the QF process includes the events of 4n emission. Such events 
account for about 1.5% of the total QF process. Generally, the QF process is regarded as 
a short process that does not have enough time to relax the degree of mass asymmetry. 
It may be unreasonable to consider that the QF process has sufficient time to evaporate 
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4 neutrons. 

It is found that among the QF process some trajectories associated with 4 neutrons 
emission enter the region of small z and their travelling times seem longer than that of 
the typical sample trajectory presented in our previous paper j22]- These trajectories 
pass the scission points where the mass number is greater than + 20 (a = 0.136) 
and less than — 20 (a = —0.136), which is classified into the QF process by the 
standard classification. In particular, it is found that the mass asymmetry of such 
trajectory changes from a = 0.67 to —0.3 in the region of small z. This QF process 
should be distinguished from the usual QF process. 

We usually consider the trajectory of the DQF process to enter the small z region, 
but due to the strong Coulomb repulsion, it goes to the large-deformation region and 
returns to the fission area without passing the fusion box. Thus, considering this 
scenario, we add a condition to the standard classification discussed in the previous 
subsection. 

Therefore, we define a new rule that a trajectory which enters the region of z < 0.4, 
where the nuclear shape resembles mono-nucleus, is to be reclassified into the DQF 
process, even though the mass fragments lie in the asymmetric region 
and A > + 20). That is to say, the DQF process is redefined as the process that 
the trajectory enters the region of 2; < 0.4 without passing the fusion box. Actually, the 
distribution of the travelling time of such a trajectory occupies the region of longer ttrav 
of the QF process. We call this classification "the classification of the fission fragment 
mass with an additional condition on z" . 

Here, we should consider the process that the trajectory goes to mass symmetric 
fission region, where the mass number is greater than — 20 and less than + 20, 
but the trajectory does not enter the region of 2; < 0.4. We name this process the 
mass symmetric quasi- fission process (MSQF). In this system, the trajectory of the 
MSQF process occupies 1.21 % of the all trajectories, and 4.55% of the mass symmetric 
fission events. We can say that the almost mass symmetric fission events come from the 
trajectory which enters the region of z < 0.4. 

By this classification, the component of the last QF process, which has a long 
travelling time, is classified into the DQF process. Such a trajectory accounts for about 
32% of the total QF process. Figure 9 shows the each neutron multiplicity distribution 
by the new classification, where each line has the same meaning as that in Fig. 3. We 
can see that the peak of the DQF process at z/„ = 2 is higher and the number of events 
increases compared with that in Fig. 3. Figure 10 shows the distribution of tfrav for each 
process using this new classification. The the travelling time of the DQF process and the 
QF process separates more clearly than that shown in Fig. 5. After ttrav = 1-2 x 10~^° 
sec, the DQF process is dominant. 

The mass distribution of fission fragments is shown in Fig. 11. The shading and lines 
have the same meanings as those in Fig. 4. We can see that the fission fragments from 
the DQF process exist in the region of mass number greater than + 20 (A = 166) 
and less than ^ - 20 {A = 126). In the regions 166 < A < 200 and 90 < A < 126, 
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Figure 9. Pre-scission neutron multiplicity given by our theoretical calculation in the 

reaction '*^Ca+^^''Pu at E*=50 McV for I = using the method of classification by 
fission fragment mass adding a condition on z. The lines denote the same items as 
those in Fig. 3. 




Figure 10. The distribution of the travelling time ttrav of each path in the reaction 
'*®Ca+^''^Pu at E* = 50 McV for / = 0, using the method of classification by fission 
fragment mass adding a condition on z. 
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Figure 11. Mass distribution of fission fragments in the reaction ^^Ca+^^^'Pu at 
E* = 50 MeV for / = using the method of classification by fission fragment mass 
adding a condition on z. The shading and hues denote the same items as those in 
Fig. 4. 

the DQF process and the QF process coexist. This means that the trajectories with 
long and short travelhng times intermingle. That is to say, in this region, the fission 
fragment has two possible origins, one is the short process wherein the colliding partner 
goes quickly to the fission area, and the other is the long process in which the trajectory 
approaches the compact shape and takes a long time to arrive at the fission area. 

3.2.4- Incident energy dependence in neutron multiplicity analysis As practice, in order 
to identify the details of the whole reaction process, we introduce the classification on the 
basis of trajectory's behavior by the theoretical calculation. When the experimental pre- 
scission neutron multiplicity is analyzed in combination with the trajectory calculation, 
the information on the neutron multiplicity can be put to practical use to clarify the 
reaction mechanism, that is to say we can distinguish the QF, DQF and FF processes. 

At a high incident energy, the situation is different. When the incident energy 
increases, the distribution of the neutron multiplicity is spread well. The neutron 
multiplicity can become a good measure of each process 1^ I28j. 

Figure 12 shows the distribution of the neutron multiplicity in the reaction 
4^Ca+24^Pu at E* = 80 MeV and E* = 160 MeV for / = in (a) and (b), respectively. 
Here, each line has the same meaning as that in Fig. 3. We use the classification 
of the fission fragment mass with an additional condition on z. With increasing the 
incident energy, the events of the QF process decrease because the mass symmetric 
fission events are dominant. The DQF and FF processes are well distinguished by the 
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neutron multiplicity at E* — 160 MeV. 

In the high incident energy case, f™""^ increases, which are about 8 and 14 at 
E* = 80 and 160 MeV, respectively. With increasing z/^"^, it becomes possible to 
identify the reaction processes using the distribution of neutron multiplicity, because 
the neutron multiplicity spreads over the wide region and indicates the characteristic 
property of each process, which is not overlapped by the fluctuations. 




Figure 12. Pre-scission neutron multiplicity given by our theoretical calculation in 
the reaction ^SCa+^'^'^Pu at (a) E*=80 MeV and (b) E*=160 MeV for I = 0, using 
the method of classification by fission fragment mass adding a condition on z. The 
lines denote the same items as those in Fig. 3 

Table 1 shows the average travelling time (ttrav) in the reaction ^^Ca+^^^Pu for 
/ = 0, using the classification method (i) the fission fragment mass (method (i)) and 
(ii) the fission fragment mass with additional condition on z (method (ii)). At E* = 50 
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MeV, in method (i), the (ttrav) of the DQF process is about 1.5 times as long as that of 
the QF process. However, in method (ii), it is about twice as long. This is because the 
DQF process in method (ii) includes the long QF process in method (i). The (ttrav) of 
the DQF process in method (ii) does not change as much as that in method (i), because 
the long tail of the DQF process mainly influences the average time. 

The (ttrav) af E* — 80 MeV and 160 MeV are shown in Table 1, in the classification 
method (ii). We can see the incident energy dependence of (ttrav)- With increasing the 
incident energy, the initial speed of the trajectory increases. After the kinetic energy 
dissipates, also the moving speed of trajectory by the fiuctuation is high. As a result, 
(ttrav) becomes shorter. The ratio among the (ttrav) of QF, DQF and FF seems to be 
approximately constant for each incident energy. 



Table 1. The average travelling time for each classification, in the reaction 

*^Ca+^''^Pu for I = 0. The methods of classification are based on (i) the fission 
fragment mass (method (i)) and (ii) the fission fragment mass with a condition on z 
(method (ii)). 





48Ca+244pu 




method (i) 


method (ii) 




E* = 50 MeV 


E* = 50 MeV 


E* = 80 MeV 


E* = 160 MeV 




(xl0-2"scc) 


(xl0-20sec) 


(xl0-20sec) 


(xlO^^Oggp) 


QF 


1.11 


0.88 


0.85 


0.73 


DQF 


1.73 


1.69 


1.58 


1.38 


FF 


4.83 


4.83 


4.00 


3.71 



4. Summary 

The fusion-fission process in the super heavy-mass region was studied on the basis of 
fiuctuation-dissipation dynamics. In order to classify the reaction process in greater 
detail, we analyzed the pre-scission neutron multiplicity connecting with the mass 
distribution of fission fragments. The neutron multiplicity depends on the travelling 
time of the trajectory from the point of contact to the scission point. We introduced the 
effect of pre-scission neutron emission into the three-dimensional Langevin calculation, 
that is to say, we combined the Langevin code with the statistical code. We applied our 
model to the investigation of the whole reaction process in the reactions ^^Ca-|-^°^Pb and 
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^^Ca+^^'^Pu at the incident energy corresponding to the excitation energy of a compound 
nucleus, E* = 50 MeV. 

In the reaction ^^Ca+^°'^Pb at E* = 50 MeV, the reaction process is simple. Most of 
the process comprises the FF process. In this case, the pre-scission neutron multiplicity 
distribution shows a single peak at = 2. 

In the reaction ^^Ca+^^'^Pu at E* = 50 MeV, the pre-scission neutron multiplicity 
distribution of the DQF process overlaps significantly with the distribution of the QF 
process, when we use the standard classification using only the information of mass 
distribution of fission fragments. Consequently, we could not distinguish between the QF 
and the DQF processes on the basis of the pre-scission neutron multiplicity. However, if 
we introduce the additional condition for the classification concerning with the travelling 
time of trajectories, we could show that the situation becomes to be resolved. The 
additional condition is that the trajectory entering into the area of z < 0.4 but not into 
the fusion box is classified to DQF process, even if the mass fragments distribute outside 
the region between — 20 and + 20. Using this classification, it was shown that 
the pre-scission neutron multiplicity distribution of the DQF process is clearly separated 
from that of the QF process. 

We also discussed how the characteristic structure of the neutron multiplicity 
depends on the incident energy, and presented that we can get more and more clear 
structure corresponding to each reaction process when the incident energy increases. 

In the next study, we would like to investigate parameters more precisely, for 
example, the friction tensor, the level density parameter and the neutron binding energy 
near di-nucleus configuration etc. We try to use the friction tensor which is derived from 
liner response theory at a low excitation energy |29 t l3( H ISTj. 
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